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Abstract 

Because of colour confinement, the physical vacuum forms an event horizon for quarks and gluons; this 
can be crossed only by quantum tunneling, i.e., through the QCD counterpart of Hawking radiation by black 
holes. Since such radiation cannot transmit information to the outside, it must be thermal, of a temperature 
determined by the strong force at the confinement surface, and it must maintain colour neutrality. The 
resulting process provides a common mechanism for thermal hadron production in high energy interactions, 
from e + e~ annihilation to heavy ion collisions. The analogy with black-hole event horizon suggests a 
dependence of the hadronization temperature on the baryon density. 

PACS numbers: 25.75.-q 
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INTRODUCTION 

Over the years, hadron production studies in a variety of high energy collision experiments have 
shown a remarkably universal feature. Prom e + e~ annihilation to p — p and p — p interactions and 
further to collisions of heavy nuclei, covering an energy range from a few GeV up to the TeV range, 
the production pattern always shows striking thermal aspects, connected to an apparently quite 



What is the origin of this thermal behaviour? While high energy heavy ion collisions involve 
large numbers of incident partons and thus could allow invoking some "thermalisation" scheme 
through rescattering, in e + e~ annihilation the predominant initial state is one energetic qq pair, and 
the number of hadronic secondaries per unit rapidity is too small to consider statistical averages. 

A further piece in this puzzle is the observation that the value of the temperature determined in 
the mentioned collision studies is quite similar to the confinement /deconfinement transition tem- 
perature found in lattice studies of strong interaction thermodynamics [3|. While hadronization 
in high energy collisions deals with a dynamical situation, the energy loss of fast colour charges 
"traversing" the physical vacuum, lattice QCD addresses the equilibrium thermodynamics of un- 
bound vs. bound colour charges. Why should the resulting critical temperatures be similar or even 
identical? 

In ref.(4), which is summarized in this contribution ( see also ref.(5)), these hadronization 
phenomena are considered as the QCD counterpart of the Hawking radiation emitted by black 
holes (BH) 0| . BHs provide a gravitational form of confinement that was quite soon compared 
to that of colour confinement in QCD 0, @] , where coloured constituents are confined to "white 
holes" (colourless from the outside, but coloured inside) . 

The main results in ref.(4) are: 

• Colour confinement and the instability of the physical vacuum under pair production form 
an event horizon for quarks, allowing a transition only through quantum tunnelling; this 
leads to thermal radiation of a temperature Th determined by the string tension. 

• Hadron production in high energy collisions occurs through a succession of such tunnelling 
processes. The resulting cascade is a realization of the same partition process which leads 
to a limiting temperature in the statistical bootstrap and dual resonance models. 

• In kinetic thermalization, the initial state information is successively lost through collisions, 
converging to a time- independent equilibrium state. In contrast, the stochastic QCD Hawk- 
ing radiation is "born in equilibrium", since quantum tunnelling a priori does not allow 
information transfer. 

• The temperature Th of QCD Hawking radiation depends only on the baryon number and 
the angular momentum of the deconfined system. The former provides the dependence of 
Th on the baryochemical potential fi, while the angular momentum pattern of the radiation 
allows a centrality-dependence of Th and elliptic flow. In particular the \i dependence of 
Tjj, will be discussed in more details in sec. 4 

THERMAL PRODUCTION PATTERN 

Let us first summarize the thermal production pattern in elementary collisions, e + e~ ,pp,pp, ■■ 
and in nucleus-nucleus scattering. 
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The partition function of an ideal resonance gas is given by 



lnZ(T) = Vj2j^Hm,T) (1) 

where di is the degeneracy factor and 4>(rrii,T) is the Boltzmann factor 

<P(m u T) = ^mfTK 2 (mi/T) (2) 
Therefore the relative abundances of the species i and j turns out 

Ni _ di4>(mi,T) 



Nj dj(j)(mj,T) 

and for transverse energy larger than T 

dN 1 



(3) 



exp--Jmf+p 2 T (4) 



dpfj. 

In elementary collisions the statistical hadronization model fits the data on the species 
abundances by two parameters: T and 7 S , that describes the strangness suppression. 

For LEP data at ^ = 91.2 Gev 0, T = 170 ± 3 ± 6 Mev and 7 S = 0.691 ± 0.053 where the 
systematic error is obtained by varying the resonansce gas scheme and the contributing resonances. 
The PEP-PETRA data at differente energies , 14 < y/s < 45 can be fitted with an average 
temperature T = 165 ± 6 Mev and -y s ~ 0.7 ± 0.05. The pp SPS data at energies y/s = 19, 23.8, 26 
Gev give T = 162 A ±1.6 Mev and 7 S ~ 0. ± 0.036 [2]. The other data for K + p and n + p scattering 
at energies close to SPS one and for pp at larger energy can be fitted by similar values. 

The fitted values of the temperature are depicted in fig. 1 [n| ■ 

Therefore there is an universal hadronization temperatute Tjj = 170 ±10 — 20 Mev which is 
independent on the species, on y/s and on the incident configuration. Moreover, also the transverse 
momentum spectra in elementary collisions can be fitted by the same value Th 0. 

In heavy ion collisions there is a new parameter which describes the finite baryon density, i.e. 
the baryon chemical potential fiB- The fits of the species abundances at high energy ( peak SPS 
and RHIC) give : 

- Tjj = 168 ± 2.4 ± 10 Mev; /x 6 = 266 ± 5 ± 30 Mev at s/s = 17 Gev for (SPS) Pb-Pb, 

- T H = 168 ± 7 Mev; fj, b = 38 ± 11 ± 5 Mev at y/a = 130 Gev for (RICH) Au-Au at y=0, 

- T H = 161 ± 2 Mev; fj, b = 20 ± 4 Mev at x/s = 200 Gev for (RICH) Au-Au. 

In conclusion, hadron abundances in all high energy collisions ( e + e~ annihilation, hadron- 
hadron and heavy ion collisions) are those of an ideal resonance gas at universal temperature 
T H ^ 170 ±10-20 Mev. 



EVENT HORIZON AND HADRONIZATION 

The idea that a thermal medium, with a kinetic thermalization by multiple partonic interaction, 
has been produced in the collisions could explain the previous phenomenon for a nucleus-nucleus 
scattering but does not work for e + e~ and hadron-hadron scattering. 

One has to look for an universal, "non-kinetic" thermalization mechanism. 

Indeed, in gravitation there is a well known example: the BH Hawking radiation has a thermal 
radiation spectrum due to tunnelling through the event horizon and the Hawking temperature is 
given by Tbh = l/8irGM, where M is the (Schwarzschild) BH mass and G is the Newton constant 
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FIG. 1: Hadronization temperature at large energy for various elementary collisions 



Therefore the conjecture [4| is that colour confinement and hadronization mechanism are in 
strong analogy with BH physics and event horizon. 

There are many reasons to believe that color confinement can be described by a color horizon 
in QCD because the theory is non linear and therefore it has an effective curved geometry 11, ]j| 
However, since we discuss the hadronization mechanism, is better to consider the Unruh effect. 

As shown by Unruh [ljj, systems with uniform acceleration, a, have an event horizon and see a 
thermal bath with temperature Tjy = ajlix. For a particle of mass,m, in uniform acceleration the 
equation of motion is solved by the parametric form 



x 



— cosh ar 
a 



— sinh ar, 

a 



(5) 



where a 
time, with dr 



F/m denotes the acceleration in the instantaneous rest frame of m, and r the proper 



v 2 dt. The resulting world line is shown in fig. 2 with the event horizon 
beyond which m classically cannot pass. The only signal the observer can detect as consequence 
of the passage of m is thermal quantum radiation of temperature 



(6) 

1/(4 GM), 



In the case of gravity, a is the "surface gravity" (i.e. the acceleration at the horizon), a 
and hence one recovers the Hawking temperature. 

In summary, the acceleration leads to a classical turning point and hence to an event horizon, 
which can be surpassed only by quantum tunnelling and at the expense of complete information 
loss, leading to thermal radiation as the only allowed signal. 
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FIG. 2: Hyperbolic motion 



On the other hand, at quantum level, it is well known that in a strong field the vacuum is unstable 
against pair production [14]]. For example, in e + e~ annihilation a qq pair is initially produced and 
when the linear potential is such that ax > oxq = 2m the string connecting qq breaks and the 
color neutralization induces an effective quantum event horizon (see figs 3,4) 





FIG. 4: Hadronization as a tunnelling process 



The qq flux tube has a thickness given by 



and the qiqi is produced at rest in e + e cms but with a transverse momentum 

1 In a 



(7) 



(8) 



The acceleration (or deceleration) associated with the string breaking and color neutralization 
mechanism turns out to be 0] 



2k T ~ \/2^a~ 1.1 GeV, 



(9) 
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which leads to 

T q = ~ ~ 180 MeV (10) 

27T V 27T 

for the hadronic Unruh temperature. It governs the momentum distribution and the relative species 
abundances of the emitted hadrons. 

Notice that the previous hadronization mechanism can be described by saying that q\ reaches 
the qiqi event horizon and tunnels to become q~2- The emission of hadron q±q2 can be considered 
as Hawking radiation. 



VACUUM PRESSURE AND BARYON DENSITY 



It is interesting to consider the extension of the previous mechanism in the case of systems with 
a net baryon number, i.e with a new "charge" which can modify the tunnelling process and the 
Hawking-Unruh hadronization temperature. 

In the BH case the effect of a total charge Q changes the Hawking temperature according to 
the formula ( see for example [l5| ) 

Note that with increasing charge, the Coulomb repulsion weakens the gravitational field at the 
event horizon and hence decreases the temperature of the corresponding quantum excitations. 
As Q 2 —* GM 2 , the gravitational force is fully compensated. The crucial quantity here is the 
ratio Q 2 /GM 2 of the overall Coulomb energy, Q 2 /R, to the overall gravitational energy, GM 2 /R. 
Equivalently, Q 2 /GM 2 = Pq/Pq measures the ratio of inward gravitational pressure Pq at the 
event horizon to the repulsive outward Coulomb pressure Pq. 

In QCD, we have a "white" hole containing coloured quarks, confined by chromodynamic forces 
or, equivalently, by the pressure B of the physical vacuum. If the system has a non-vanishing 
overall baryon number, there will be a Fermi repulsion between the corresponding quarks, and this 
repulsion will provide a pressure P (/i) acting against B, with u denoting the corresponding quark 
baryochemical potential. We thus expect a similar reduction of the hadronization temperature as 
function of \i. To quantify this aspect let us consider, at T = 0, the Fermi pressure, P = dgfi 4 /2Air 2 , 
where d g is the degeneracy factor, versus the QCD vacuum pressure due to the gluon condensate 
B =< (a s /i{)G 2 > which is a decreasing ( largely unknown) function of the baryon density. By 
following the analysis of ref. [lfj for the dependence of the gluon condensate on the baryon number, 
the critical density, where B balances the Fermi pressure, turns out about n c = 5 — 6no, with no 
the nuclear saturation density. 



A slightly different result is obtained in the description of deconfinement by percolation 17( : 
n c ~ 4n . 

In the intermediate region, where both T and \x are finite, we want to compare the effect of 
the Fermi repulsion to the vacuum pressure through the Hawking-Unruh form, i.e., we replace 
Q 2 /GM 2 in eq.(ll) by O/W)) 4 , giving 



T ^ )/To ~ (i + vi-w^ - (12) 

The resulting behaviour of T(/z) is shown in Fig. 5 for hq corresponding to the previous n c ~ 4 — 6 
riQ. In the same figure is shown the results obtained by using (h/hq) 2 rather than (fi/fio) 4 in 
eq.(12). In both cases the function remains rather flat up to large value of /i. 
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FIG. 5: T dependence on jj, by eq.(12) with (/j,//j,q)" for n = 2,4 



Clearly this approach is overly simplistic, since it reduces the effect of the additional quarks to 
only their Fermi repulsion. A more general way of addressing the problem would be to introduce 
an effective /i-dependence of the string tension. 
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